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Abstract: We study the link between neutrino oscillations and leptogenesis in the 
minimal framework assuming an SO (10) see-saw mechanism with 3 families. Dirac 
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1. Introduction 

The evidence for a total asymmetry between baryon and anti-baryon densities in our 
surroundings is an indisputable fact of life. This could easily be accounted for by 
assuming that the initial conditions for our universe are such that this asymmetry is 
of order one. However, there is every reason to believe that the thermal history of the 
universe can be traced back in time using known particle physics up to temperatures 
of maybe 100 GeV, certainly high enough to massively produce quark-antiquark 
pairs. At these temperatures, big-bang nucleosynthesis requires Q that the adiabatic 
invariant ratio 

ys = n £ -n 1 = 0710 _ 10 

s 

meaning that there be about one extra quark for 10 10 quark-antiquark pairs. Avoid- 
ing such extreme fine-tuning by a dynamical mechanism able to produce this number 
out of an initially symmetric configuration is the purpose of baryogenesis. 

Many mechanisms of baryogenesis have been proposed in the past 30 years (see 
e.g. || for a review including the most ingenious and exotic ones), but most rely 
on ad-hoc new physics. We certainly would prefer to explain the baryon asymmetry 
from solid, experimentally tested physics. The standard model of particle physics 
satisfy this criterion at a desperately high level of precision, and it was shown || 
(see [|J for a review) to satisfy in principle all the Sakharov|J conditions required 
to produce a baryon asymmetry. However, what looked like a very small number 
for initial conditions, now appears too large for the pure standard model to achieve. 
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At least, CP violation beyond the known CKM phase must be added to resist the 
strong GIM cancellations in the hot plasma 



Furthermore, the only possibility to change baryon number in the standard model 
is via sphaleron processes which should freeze abruptly below the electroweak phase 
transition to leave a net asymmetry 0. Since the transition gets weaker when rising 
the mass of a single scalar doublet, extra non-standard scalars are invoked to counter 
the rising of the experimental lower bound on the lightest scalar. Supersymmetry 
naturally provides a lot of well-motivated extra scalars, but getting a strong enough 
transition with inn ~ lOOGeV requires a huge gap between the right-handed stop and 
all other sfermions, which may seem contrived, and even so, could not protect E-W 
baryogenesis against a « HOGeV Higgs bound ||. Unless a scalar is soon discovered, 
we thus seem to be lead back into the pre- [[[(]] situation, where baryogenesis was a 
footprint from extreme high-energy physics, with little chance of an experimental 
cross-check. 



However, another experimental signal for non-standard physics has since then 
developed into a quantitative and solid field, namely the evidence for neutrino oscil- 



lations which strongly points toward small but non-zero neutrino masses [TT|]. If these 
masses are of the Dirac type, right-handed neutrinos must be added, but we are left 
with the puzzle: why are neutrinos 10 10 times lighter than charged leptons? In our 
mind, the smallest theoretical price to pay for resolving this puzzle is to keep the 
right-handed neutrinos, but give them large Majorana masses ~ 10 10 ~ 16 GeV: this 



is the celebrated see-saw mechanism ||12||,||13||, ||14|,fl5|. Such a high Majorana mass 
breaks lepton number at a slow enough rate, and in a very indirect way, one can say 
that neutrino oscillations provide a leptonic analogue of the baryon number violation 
looked for unsuccessfully in proton decay searches during the last decades. 

With this theoretical prejudice, we are thus lead to take the existence of heavy 
right-handed neutrinos for granted. This makes leptogenesis [pl| an extremely natu- 
ral mechanism to consider for producing the baryon asymmetry. Indeed, the decay of 
right-handed neutrinos can easily leave a CP-odd lepton asymmetry, which standard 
sphaleron processes can convert into a baryon asymmetry. It is then most interest- 
ing to see how closely can this asymmetry be related to tested or testable neutrino 
oscillation physics, and how much asymmetry can be produced. 

To study these questions, we first review in section § the generic leptogenesis 
mechanism, to fix notations and give the formula for the asymmetry in terms of 
neutrino masses. Then we focus in section |ona generic class of 3 families 5*0(10) 
models offering a computable connection between leptogenesis and neutrino oscilla- 
tions. In section |], we analytically work out this connection in a toy see-saw model 
with only 2 families, to help and understand the results presented in section |^. 
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2. Leptogenesis 



One of the most attractive scenario to explain the baryon asymmetry of the universe 
is leptogenesis. In such a scenario, a primordial lepton asymmetry is generated by 
the out-of-equilibrium decay of heavy right-handed Majorana neutrinos N{. The 
possible excess in lepton number L is converted into a baryon asymmetry B by the 
sphaleron processes. From relations between chemical potentials this conversion is 
given by (see |I?|,|I8| and references therein): 



SN F + i N H 

22N F + 13N H K ' v ; 

where Np and Nh are respectively the number of families and Higgs doublets. 

At low energies, right-handed neutrinos (iVj) are gauge-singlets, and only appear 
in the Lagrangian via their Majorana masses (Mj) and Dirac Yukawa couplings (hu) 
to the Higgs doublet (H) and lepton doublet (Li) 

C = £ M + C Y = ^{NiT c M l N l + L~ l i(j2H*h h N i + h.c (2.2) 

in the right-handed neutrinos eigenbasis (0 < M% < M 2 < M 3 ). By the see-saw 
mechanism, light neutrinos then get a Majorana mass-matrix : 

m v = v 2 /i.Diag(— , — , — ).h T (2.3) 
BK Mi M 2 M 3 J V ; 

The 3 Sakharov conditions are easily met by the iV-decays. Indeed, the right Ma- 
jorana mass fulfils the first condition by violating lepton number L, manifest in the 
simultaneous existence of both decays channels iVj — > H + I and Ni — > W + I. 

Second, the CP violation condition can be met by the interference between tree- 
level and e' (vertex) or e (wave function) one-loop corrections. Indeed, CP asymmetry 



in the decay (in vacuum) of the right-handed neutrino iVj is [19, 20|: 



5 _ E ; r(iv t ->/ + g)-r(iv t ->I + gt ) _ 



with 



. 1 MEW MjMj , 1 ImEW f( M[> 

H ' j Svr Y.Ahu\ 2 Mj-M?' ^ Svr E J^| 2 Mf' 



where f(x) = y/x [l — (1 + x) ln(i^)] . Defining 

Aij = (rfh) tJ (2.6) 
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the virtual Nj contribution to the Ni decay asymmetry becomes for a large hierarchy 1 
between heavy Major ana neutrinos Mi <C Mj 



1 Im(4) M t , 1 Im(4) M, 

"' 8tt An Mj ' ^ ~ 16tt An \\ , 1 ' 



while in the opposite limit Mj ^> Mf 

1 M4) Mj , 1 Im(A|,) M 3 | M 3 
Stt Aii M' JJ ~4vr A„ Mi { Mi } - { ' 

Concerning the third Sakharov condition, let us first assume Ni to be thermally 
distributed in the early universe when T > Mj. While cooling, the equilibrium 
distribution starts feeling a strong Boltzmann suppression when T « Mj. If the 
vacuum decay rate = ^-AuMi was much slower than the expansion rate at that 
time H(Mi) = each N would decay essentially in vacuum and starting 

from 21 « 1 , we would get |23) Y* = Lj/s ~ ^/t?*. However, decay rates are not 
necessary that slow, and if the universe is still hot by the time Ni decays, inverse 
decays dilute Si roughly into 5i/K~i with : 

_ r JVi (T = M) „ 1 , 

Jvj — „ ,„ : , ; ~ ~ — : — ~ ^=/iii" 



#(T = M) 8?rl.66 M 

characterising the efficiency of (inverse-) decay processes when N becomes non-relativistic. 
It is convenient |24] to rewrite this dimensionless parameter in terms of an effective 
see-saw like mass 

- . v 2 A ti 
m; = —-r- = Kim 

Mi 

with a critical value m* = a/512^71- 5 /90 v 2 / M p* = L08 10" 3 eV in the Standard 
Model. Below this critical value, decays happen like in vacuum, but the Yukawa 
couplings are too small for the iVj population to reach equilibrium, and the result 
depends on the choice of initial conditions. If we assume no initial population, there 
is also a suppression for rhi < m*, and the final asymmetry 

S 

Yi = —d{m h Mi) 
9 

contains a dilution factor d(m) that is only close to 1 for m ~ m*. For Mj < 10 8 GeV 
where 2 body scatterings can be neglected, the solution of the Boltzmann equations 
(fig. 6 of P3| ) can be reasonably fitted by 

1 



d{m) 



(m/10m*)-°- 8 + (m/0.28m*) L25 

which is compared with the usual approximation p3| d(K > 10) ~ 0.3(i^ln ' 6 K)~ 
on figure [I] (see also |p6| ). 



1 The expressions ( |2.5| ) anyway require resumming in the degenerate limit, as the asymmetry 
should then vanish instead of diverging pH E(X flty, E3|. 
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Figure 1: Fits of the dilution factor d computed by |2f| in the SM for M\ = 10 8 GeV 
(red dots): the plain blue curve is our formula (||), to be compared with: 1) the usual 
approxmation 0.3(-fT In 0,6 K) -1 for large K from p3], extended to lower if' s following 



(green dashed curve); 2) the expression of 28], more appropriate for larger M\ (purple 
dotted curve). 



Decay of the lightest Ni usually gives the latest and dominant contribution, 
because comparing (|2.7p and (|2~8| ) shows that other A^-'s are suppressed by A^/A 2 -- 
(typically 1CT 10 ) and further diluted by scattering processes (see for instance [^9| 
for a discussion of such processes). Using Yb ~ — 1/3 Yl from sphaleron conversion 



(2.1|), the final baryon asymmetry we use in section 5 is then: 



10 10 -3Mi x-^ Im(A?,-) 

y BW = 10 » y fl „ E „ o. 7 (2.9) 

J — 2,3 



At this level, no constraints from neutrino oscillations encoded in ( |2.3|) can be 
drawn because the Yukawa coefficients h are free parameters: for any set of pa- 
rameters (h, M) providing good oscillation physics but bad Ybw, we can rescale 
h — ► h = h/y/Ysio and M — > M' = M/Ybiq to get the same light neutrinos, with 
^sio ~ 1- Notice that K and rh are invariant under such rescaling, so that the only 
violation to this scaling law comes from the Mi dependence of the dilution factor d. 

To set the M\ scale of interest, it is instructive to make the further approximation 
d(l < K < 10 3 ) ~ 0.1/ K and assume domination of the Ni decay with a virtual A^ 2 
exchange. The result for the asymmetry is then 

0.7 x 10- 10 Im(A? 9 ) Mi M x , 
Ybio « 7= y^-T^-TK^Tr « 0.7 (2.10) 



ftf A 2 U M 2 10 10 GeV 
so that Mi « 10 10 GeV is the natural scale for leptogenesis. 
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To draw any link between neutrino oscillations and leptogenesis, we thus need 
further theoretical assumptions, of which we now describe our minimal set. 



3. SO(10) 

5O(10) D SU(3) x 517(2) x U(l) 

The iSO (10) Grand Unification group is a very natural framework for studying 



leptogenesis. Indeed |30], |3T|, it is the smallest GUT group including right-handed 
neutrinos together with usual fermions in matter multiplets : 



16 



(3, 2) + (3, 1) + (3, 1) + (1, 2) + (1, 1) + (1, 1) 



N 



Furthermore, it includes (B — L) as a gauge symmetry, whose breaking gives a natural 
explanation for right-handed neutrino Majorana masses. Finally, leptons and quarks 
masses are related which fixes the neutrino Dirac Yukawa h. For instance, if Dirac 
masses originate from a single scalar in the vector (10-dimensional) of 50(10) : 



(3.1) 



we obtain the (unphysical) low energy relation : 

nid = 3m e = m u = 3m u oc f 10 v 10 



where the factor of about 3 comes from the extra color loops affecting quarks but 
not leptons. Adding a second (or more) scalars in the same 5O(10) representation : 



Y 



(3.2) 



lifts the degeneracy m u = rrid while keeping m u k, ?>m v i.e m v rs lMeV — > 50GeV. 

To get more realistic neutrino masses, an extra scalar 126 in the 126-representation 
with v.e.v v\26 in the 5?7(5)-singlet component can be invoked to break (B — L), pro- 
viding the Majorana mass term : 



£m = f^v 12S N-N b + h.c 



(3.3) 



Thus, in the framework of 5O(10), Dirac masses are naturally determined by up- 
quarks masses. 

Aij expression 

Writing the relevant Lagrangian term for the leptons we get after SU (2) breaking : 



C 



L a x 



fab^Rb + 



L a X 



f: b N b + \mv l2& f™N b + h.c (3.4) 
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so that M a b = /^ 6 fi26 with a and b being components in some arbitrary family 
basis. Decomposing the Yukawa couplings into mixing and real masses eigenvalues, 
we write : 

f d = ^Ul d .Di&g(m d ,m s ,m b ).U Rd (3.5) 

/" = ^:t / L-Diag(m u , m c ,m t ).U Ru (3.6) 
ov 

f 126 = ^^/•Diag(M 1 <M 2 < M 3 ).Ul (3.7) 

The previous hu are up- Yukawa matrix elements between charged leptons mass (and 
flavour) eigenstates (/;) and right-handed neutrinos mass eigenstates (Ni). Therefore: 



so that 



U Ld .f u .U M (3.8) 

—U Ld .Ul Dia,g(m u , m c , m t ).U Ru .U M (3.9) 
6v u 

^v£ KM .Dia,g(m u ,m c ,m t ).V R (3.10) 

A H = Q^ v l D M™l,rn 2 c ,rr$).V R (3.11) 

only depend on the unknown V R mixing. But the light neutrinos mass matrix is 
constrained by neutrinos oscillation physics: 

m v = ADiagf — , — , — ).h T (3.12) 

= C/.Diag(mi, m 2 , m 3 ).U T (3.13) 

The unknowns (V R , M 12 ,3) can thus be determined from measurable quantities (U, rni^s) 
by diagonalizing the matrix : 

M* R = -Diag(m u , m c , m t ).[/ e// .Diag( — , — , — ).Uj ff .Diag(m u , m c , m t )(3.14) 
9 mi m 2 1713 3 1 

= ^.Diag(M 1; M 2 , M 3 ).V% (3.15) 

where 

U eff = V CKM .U (3.16) 

Parameterisation & experimental inputs: 

Using the standard parameterisation ||32|| for V : 

C12C13 S12C13 s* 3 

V CKM = V23.V13.V12 = ( -S12C23 - C12S23S13 C12C23 - S12S23S13 s 23 ci3 I (3.17) 

S12S23 - C12C23S13 -C12S23 - S12C23S13 C23C13 , 
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we use the following values for the quark mixing angles: 

s g KM = 0.041 

sg KM = 0.0036 e iS ° KM ; 8 CKM * 
sg KM = 0.223 



As for light neutrinos ||33||, we label mass eigenstates so that the most constraining 
experiments each fix a definite mixing angle and use a similar decomposition with 
three complex mixing angles, and Majorana phases: 

U = U 23 . Uu . U 12 .Diag(l, e^ 1 , e^ 2 ) (3.18) 

atm Chooz sun 

s l2Catm C sun S2^Si2, C sun C a tm -5l2^23^13 ^23 ( -'e3 
S12S23 — C sun C a tmSl3 —C sun S23 ~ S 12 C f m Si3 C a t m C e 3 

where the complex parameters are: 

a — a e^ 13 





(I \ 




e^ 1 




V0 e^ 2 / 



S23 — s a t m e 



Squared mass differences are then fixed (up to a sign) by the oscillation parameters 

Am2 sun < Am ltm- 

m\-m\ = Am 2 sun (3.19) 
m\-m\^m\-m\ = ±Am 2 atm (3.20) 

In both cases, is singled out by atmospheric oscillations, but according to the 
sign, it can either be heavier 

mi < m 2 <C m3 (standard hierarchy), 
or lighter than the others 

m 3 <C rrii < m 2 (inverted hierarchy). 

We take the lightest mass (mi or TO3) as the free parameter that cannot be fixed by 
oscillations. For values of this parameter larger than Am atm , we get two different 
realizations of the degenerate case. In the following, it will vary between 10 _7 eV 
and leV: going below this range requires huge right-handed masses, typically above 
GUT or even Planck scale; going above does not help to increase the asymmetry, and 
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generally violates neutrino- less double beta decay experiments |34], ^] (see however 
[p6|| for a study of possible cancellations using Majorana phases). 

Other free parameters are the six phases: S (from V CKM ), 5u, 613, 623 (from 
the complex mixing angles 2 in U), and the two Majorana phases 01,2. A careful count 
of the 6 phases in see-saw models can be found in reference [[37] . 



With the above parametrisation, the experimental constraints from oscillations 
lead us to take the following values in the subsequent numerical analysis: 



Atmospheric neutrinos [32]: Am„ im ^310 eV with nearly maximal mixing: 



sin 2 20 atm = 0.92 1 so that s atm = 0.6 -> 0.7. 

CHOOZ reactor experiment [|38|] : for the above value of Am 2 tm , |C/ e3 | = s e s < 
0.18 at 95% C.L. 



Solar neutrinos |39|]: we will consider 4 possible resolutions of the solar neutrino 
deficit: 

1. large mixing ang le MSW oscillations (LMA) with Am 2 sun « 5 10~ 5 eV 2 and 
tan 2 snn = 0.42 G [0.15,0.9] at 95% C.L. corresponding to s sun = 0.54 G 
[0.36,0.69] ; 

2. small mixing ang le MSW oscillations (SMA) with Am 2 sun « 5 10" 6 eV 2 
and t&n 2 8 sun = 1.5 10~ 3 G [6 10~ 4 , 3 10~ 3 ] at 95% C.L. corresponding to 
s sun = 0.039 G [0.025,0.055] ; 

3. LOW with Am 2 un w 7.9 10 _8 eV 2 and tan 2 8 sun = 0.61 corresponding to 

Ssun 0.62 , 



4. vacuum oscillations (VAC) with Am 2 Mn ps 4.6 10 10 eV 2 and tan 2 6. 



sun 



1.8 so that s sun = 0.8; 
Although the LMA solution is favoured by the SNO data HD[ and confirmed 



by the latest KamLAND results 0, we still consider the SMA, LOW and VAC 



situations for completeness and illustration. 
4. Toy see-saw (2 flavours) 

Before reviewing the results in a full case with 3 families, it helps intuition to work 
out a 2 families exercise. In particular, equation ( |2.9|) suggests that leptogenesis 
prefers right handed neutrinos with large mixing and similar masses. This is not 



easy to achieve by see-saw with large Dirac hierarchies and mixing [42, 43 . Starting 



from the see-saw formula ( |3.14| ), dropping the first family and introducing the Dirac 

2 Notice that we take 612,13,23 € [— tt/2], so that s stm , e3 , atm € [-1,1] and c su „, e3 , Qtm G [0,1] 
instead of being restricted to the first quadrant like the usual CKM angles 
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hierarchy parameter r = m c /m t ~ 10 2 for instance and the light neutrinos hierarchy 
parameter m = 7712/1713, we get in the case of real effective mixing parametrised by 

e/j 
23 



s = sin 6, 1 ■ 



9 \ m t J \ — s c 
m t ( r 2 (c 2 + s 2 m) rcs(m — 1)\ 1 




9m3 \ rcsim — 1) c m + s J m 



(4.2) 



For a large inverse hierarchy (w.r.t. Dirac hierarchy at small angle) m > 1, we 
get for small r,s < 1 (in units of m 2 /9m^) 



r 2 s 2 + r 2 /m rs\ ( c R —s R \ ( M\ \ / c R s R 
rs 1 J ~ \s R c r ) ' V mJ ' \-s R c R 



(4.3) 



so that M 2 « TtM r « 1, M a « DetM fl = r 2 /m < TrM fl while the right-handed 
mixing ~ r 2 s 2 is strongly suppressed by the Dirac hierarchy r <C 1. The ratio 
appearing in the leptogenesis ( |2.9| ) is 

A12 _ c R s R ^ _s 
A u r 2 c\ + s\ r 

and the the asymmetry ( p.lOj ) goes in this case like 

mo A 2 X M 2 m 2 

which is strongly suppressed and dominated by the smallest values of m possible. But 
in this degenerate region m ~ 1, another type of suppression arises as vanishes 
together with the off-diagonal terms in ( |4.2| ) . 

To hope for larger results, we must thus turn to the usual hierarchy limit m < 1. 
Assuming again r,s < 1, M R in ( f4.2| ) then becomes 

r 2 —rs \ 1 



— rs s 2 + 771 J 771 

By furthermore assuming Mi <C M 2 (i.e. r 2 m <C (m + r 2 + s 2 ) 2 which turns out to 
be satisfied everywhere, except when s 2 <C r 2 m m (see figure ||), we arrive at the 
useful analytic expressions 

1 

M 2 « Tr M R « — (m + r 2 + s 2 ) (4.4) 
m 

DetM^ r 2 

Mi « « (4.5) 

1 Tr Mr m + r 2 + s 2 1 J 

— TS 

-A21 = s K c B — — — (4.6) 

m + r z + s z 
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Figure 2: Right-handed masses Mr in arbitrary units, as a functions of light neutrino 
masses ratio m = 1712/1713, for 3 values of light neutrino effective mixing s = sin(#23e//) 
10 . 10~ 2 , 10~ 3 and a fixed value of the Dirac masses hierarchy r = m c /mt = 0.01. For 
s = 0.1, we also show the effect of the Majorana phase 4> m = 0.99-7T (red curve) leading to 
degenerate Majorana masses. 



From these expressions, it is clear that a cross-over occurs at values of m around the 
largest of s 2 or r 2 , but unlike the level-repulsion effect usual in quantum mechanics, 
the "levels" M 12 do not need to be close for the cross-over to happen, as can be 
seen in the s > r situation. In this case, the levels do not asymptotically tend to 
the no-mixing situation away from cross-over region. These results are most simply 
summarised in the plots |2|, drawn without approximation for r = 10~ 2 . 

Since An = r 2 c 2 R + s 2 R , the ratio Ai 2 /Au only takes a simple form for small 
mixing sj C r (large enough m or small enough s), in which case the expression 
coming into the baryon asymmetry is simply 

growing like m cx M^ 1 for small m, reaching a maximum around r 2 + s 2 and de- 
creasing like m -4 beyond. We see in figure |3] that these approximations capture 
the general shape of the asymmetry, except for extra suppressions in the regions 
m ~ 1 (left-handed degenerate ), and m < s if s 2 < r 2 {An ^> r 2 ), where the 
approximations used break down. 

Of course, the true asymmetry Y BW vanishes unless A 12 gets a CP-violating 
non-zero phase. Specialising (|3.18| ) for 2 flavours leaves 2 phases, 4> m = 2(0 2 — 4>i) 
and S23, which respectively complexify the parameters m = |m|e J ^ m , s = \s\e zS23 
and sr = \sn\e tSR in the previous discussion. It is a bit tedious to check that the 
above approximations (|4.4H4.q) for Mi <C M 2 and \A^\ can be extended to complex 
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Figure 3: Log-Log plot of the asymmetry Ygio (in arbitrary units) as a function of 
m = 1712/1713 with the same conventions as in figure |2[ For each curve there are 2 twins: 
the red spiky curve is the exact version of (4.9) computed for the phase maximising its 



value, while the smooth blue curve is the exact version of (4.7) for real positive parameters. 



parameters, and that the relevant phase is given by 

A 12 = e- i5 «\A 12 \. (4.8) 
The asymmetry ( [2.10| ) in the same approximation as ( ^4.7| ) is then actually 

y^W.^.M,^^ (4 . 9) 
A\ x M 2 \m + r z + s z \° 

which has the same qualitative features as the previous expression , with a 
further suppression by the factor | sin(— 2<5r)| < 1, and a possible enhancement when 
the denominator cancels, i.e. outside the range of validity of the approximations. 
Figure ||| shows that for a maximising value of the phase 6r, the correct asymmetry 
Y BW is close to Y^ r worked out without phases, after replacing Im(v4 2 2 ) by \A\ 2 \ = 
A\ 2 . This approximation can be a useful guide for finding upper bounds in the full 
3 flavours case, where the large number of phases makes it non-trivial. 

To summarise, the lesson learned in this 2 flavours see-saw exercise is that 

1. the right neutrinos mass ratio M\jM 2 is maximised for a ratio of light neu- 
trino masses m = m\jm2 which is either the squared ratio of Dirac masses 
r 2 = (mDi/iriD2) 2 or the squared sine of the light effective mixing angle s 2 = 
sin 2 O e ff, whichever is largest, 

2. the baryon asymmetry can be maximised for a similar or slightly larger value 
of m, 
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3. Majorana phases cancellations can lead to further enhancements around that 
same value of m (red curve on figure ||), but these can only be trusted after a 
careful resummation of this degenerate Mi ~ M 2 case; 

4. for small mixing, the "inverted hierarchy" m\ ^> m 2 gives much smaller results 
than the "standard hierarchy" m\ <C mi. 

Finally, we should also notice that this 2 flavours exercise was based on the 
approximate form ( |4.7| ) which, as ( |2.10| ), is only valid for K > 1, and overestimates 



the asymmetry for lower K. For the results below we use the relation ( |2.9| ) which 
includes more realistic Boltzmann dilution. 



5. Results with 3 flavours 

Let us now discuss the more realistic situation of 3 light neutrino flavours. It turns 
out the baryon asymmetry is generically too small to be useful. Before concluding 
that 50(10) see-saw leptogenesis is excluded, we must look for regions of parameters 
where the asymmetry can be maximised. As seen in the previous section, standard 
hierarchy is a better candidate in this respect. We will pass in review the various 
solutions to the solar neutrino deficit, starting by vacuum oscillations (VAC). 



Vacuum: Right handed masses 
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Figure 4: Masses of right-handed heavy Majorana neutrinos for vacuum oscillation of 
solar neutrinos (Am sun = 2.1510 _5 eV). The dash-dotted curves completely ignore mixing 
(U e ff = 1), the plain curves are for real CKM mixing, maximal atmospheric and solar 
mixing but vanishing s e ^ while the dashed curves tune s e 3 to maximise Mi/M%. 

In view of the previous section, a good starting point to maximise the baryon 
asymmetry is to try and make the lightest right-handed masses Mi and M 2 as close 
as possible to each other. If there were no mixing, U e ff in ( pM4j ) would be diagonal 
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and right-handed masses would follow the dash-dotted curves of figure [|. For large 
mi (degenerate limit) they draw 3 parallel lines separated by large Dirac mass hi- 
erarchies. With decreasing mi, all Mrs grow until m 1 = Am atm w 10 _1 eV, below 
which M 3 levels off at m^/9Am atm ph 10 14 GeV. Meanwhile M 12 continue growing 
together until mi = Am sun ph 10 _5 eV, where M 2 stops at about the same value. 
Further decreasing mi then only affects Mi, which becomes degenerate with M 2 for 
extremely low mi ~ iO~ n eV. 
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Figure 5: The baryon asymmetry Ysio for LOW and VAC oscillations of solar neutrinos. 
The plain curve is the vacuum case for vanishing s e 3. The gain from tuning the value 
s e 3 = —0.16 is shown for the VAC (dashed) and LOW (dash-dotted) cases, taking for each 
case the maximising value of the phase #13. 



If we now turn on the largest allowed atmospheric and solar mixing, we get the 
plain curves in figure f|. Decreasing mi from the unaffected degenerate limit, M 3 
starts again levelling off at mi = Am atm . But the maximal atmospheric mixing 
immediately induces the type of "level crossing" seen in the previous section, which 
effectively exchanges M2 and M3. M2 thus levels off at a much lower value ph 10 10 GeV 
which offers a better possibility for Mi to catch up. However, despite keeping s e3 = 0, 
CKM mixing in eq. ( |3.16| ) still induces a non-trivial |si 3e //| ~ 0.16 which stops the 
growth of Mi below mi = \siz e ff\ 2 Am atm ph 10 _3 eV. M 2 on the contrary starts 
growing again until it is hit by the solar mixing at mi ph Am sun . We thus expect a 
maximum baryon asymmetry around mi ph 10~ 3 eV in this case. A plot of the full 
3 flavour asymmetry (plain curve in figure |5|) confirms this expectation, and further 
shows that this maximum is at least 5 orders of magnitude too low. The shape of 
the asymmetry nicely fits the picture derived in the previous section, once we recall 
that for mi < Am sun both Mi and M 2 stay constant. It is also worth noticing that 
at this stage, the only source of CP violation is S CKM , and that it drives the right 
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sign for the baryon asymmetry, albeit with a small amplitude. 
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Figure 6: Sensitivity of the baryon asymmetry on the matrix element U e3 = s e3 e lSl3 most 



constrained by CHOOZ. The levels of Ybio shown are 1(T 6 , 1(T 5 
only obtained on the darkest dots around (s e3 = — 0.16, #13 = 
zoom on the right. 



...10 2 , the latter being 
0.0) at the centre of the 



To get larger results, we may use the freedom left by the CHOOZ experiment 
to play with s e3 . Indeed, for s e3 ~ —0.16 (dashed curve in figure |), si 3e // ~ so 
that Mi is decoupled from the value of m 3 , and continues rising up to m^/9Am sun ps 
10 8 GeV. This may give a rise in the asymmetry down to mi ~ 10 _5 eV. Notice 
that for such a value, K ~ 1 which minimises dilution effects, M\ ps 10 8 GeV which 
is safe w.r.t. gravitino bounds on inflation in SUSY versions of this scenario, and 
M 3 ps 10 16 GeV which is marginally compatible with the GUT scale. It is worth 
detailing where the special value s e3 ~ —0.22 s atm comes from. In the expression of 
U e ff, we may in a first approximation neglect all CKM mixing except the Cabbibo 
V\2 KM ■ Then keeping terms at most linear in 8y£ and s e3 , we may write: 



U, 



eff ~ Vg KM U 23 U 13 U 12 

= u 23 u 13 [ul 3 ul 3 vg KM u 23 u 13 
« u 23 u 13 v; 3 v; 2 u 12 

= U 23e ff Uis e ff Ui 2e ff 



u 
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(5.1) 
(5.2) 
(5.3) 
(5.4) 



with the commutator V matrices parameterised like fl3.17| ) with s' 12 ~ c atm S\ 2 and 
s' 13 ps s atm e 4<523 Si 2 KM . Cancelling Si 3e // then requires to satisfy the complex equation 



U, 



c3 



, atm t 



-1S23 Q CKM 
5 12 



-0.16e 



-i$23 



(5.5) 



which is illustrated for #23 = in figure || One sees the extreme tuning on both the 
modulus and the phase of U e3 , needed to achieve an enhancement in the four dark 
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central dots of the zoom. Even the sign of the asymmetry violently flips in the four 
quadrants of this zoom. 
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Figure 7: The baryon asymmetry Ybio f° r MSW oscillations of solar neutrinos with various 
mixing angles. The plain curves are for vanishing s e 3 with no phase except 5 . The 
dashed SMA (resp. LMA) curves are for Sis = —0.02 (resp. —0.1) . Only the particular 
dash-dotted curve (s e 3 = —0.16, s sun = —0.156, #13 = 0.005) can nearly produce enough 
asymmetry, but it can not account for the solar neutrino deficit. 



The previous discussion can be repeated for the LOW solution, but the larger 
value of Am 2 sun leaves less room for the asymmetry to grow before levelling off. As 
seen in figure [5], the maximum that can hoped for is 2 orders of magnitude too low. 

If we now turn to MSW solutions of the solar neutrino deficit, we get a still larger 
value of Am sun , and there is too little room for Mi to grow between mi ps Am atm and 
mi ps Am sun . However, it can grow for smaller mi, provided both Si3 e // = sueff = 0. 
This is in principle possible, but now requires very special values for both s e 3 ps —0.16 
(as previously) and s sun ~ — s^ M c atm ~ —0.16 (see eq. |5T3| ). This last value of s sun 
is incompatible with any possible MSW solution, unless taking away from its 

quark-lepton symmetry value, which is illustrated as a function of mi and s sun on 
figures and ||. 

With a standard mass hierarchy, we thus see that vacuum oscillations are the 
closest to account both for the solar neutrinos deficit and the baryon asymmetry of 
the universe. With the same line of reasoning, it is easy to see that with inverted 
hierarchy (m 3 <ti nil <ti m 2 ) the asymmetry is even smaller, because of a larger Mi — 
M2 gap coming from the bound Mi < m„/6Am a(m . Indeed, as the free parameter 777,3 
is lowered below Am alm , two Mr eigenvalues must now stop growing instead of one. 
There are then two possibilities: if Mi is the only one that grows in the absence of 
mixing, then maximal atmospheric mixing stops its growth once m.3 < s atm Am atm . 



- 16 - 




-1 -0.5 0.5 1 

Ssun 



Figure 8: The baryon asymmetry Ybio f° r MSW oscillations of solar neutrinos (LMA: 
lower red, SMA: upper blue), as a function of s sun . The maximum is nearly large enough 
(and sharply peaked) for mi ~ 10 -5 if s e s = —0.16, but lies just between the small and 
large mixing angle rectangles. 



If on the other hand M\ stops growing without mixing, then mixing won't help, as 
it can only induce "levels repulsion", not attraction. In both cases, M\ is at most 
10 4 GeV, and Mi/M 2 cannot exceed m^/ra^ ~ 10~ 6 which makes the asymmetry too 
small. 

Let us now review the effect of the various phases for trying to reach a sufficient 
asymmetry. In the discussion above, we have only switched on the minimal phases 
necessary to get a non-trivial result, namely S CKM which is unavoidable in our quark- 
lepton symmetry set-up, and 813 which was necessary to enforce the cancellation ( |5.5| ). 
From this equation, it is clear that if #23 were non zero, it would merely shift the 
optimal value of #13, the asymmetry being a function of the combination 5 13 — 5 2 3 
in the region of the low m x peak. An non-zero 612 would induce a much smaller 
shift, through the sub- leading corrections to ( |5.5|) . Less trivial are the effects of the 
Majorana phases 0i ; 2- As suggested by equation ( [Of ), these can lead to cancellations 
and near degenerate right-handed masses for a fine tuned value of one of the phases 
and of m\. However, the peak in the decay asymmetries ei 2 and ei 3 tend to cancel 
each other, and get further reduced by a similar peak in K%, so that no enhancement 
of the final asymmetry can be produced for the large mi « 10~ 3 eV where it could 
occur. 
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6. Conclusions and perspectives 



In this work, we have studied the possible relations between neutrino oscillations 
and the baryon asymmetry of the universe produced by equilibrium decays of right 
handed neutrinos, assuming an SO (10) see-saw mechanism for neutrino masses. As 
explained, this is a minimal predictive set of assumptions to study such relations. We 
find that the produced asymmetry is generically six orders of magnitude too small 
because of the huge Dirac masses hierarchy assumed (r 2 = m 2 /m 2 ~ 10 -6 ). The 
closest one can come to the observed asymmetry in this framework is at least one or- 
der of magnitude too small and requires vacuum oscillations of solar neutrinos, with 
standard hierarchy and very specific values for the least constrained neutrino param- 
eters: 1) the lightest mass (mi ~ Am sun « 10 _5 eV), and 2) the heavy component of 
the electron neutrino (|s e3 | = |s atm |.| sm(9 Ca bbibo)\ ~ 0.16). 

This result calls for several extensions and refinements. It would first of all be 
more satisfactory to get the right-handed mass relations we tend to need out of some 
theoretical mass model for all particles. This seems quite challenging at the moment. 
Nice family symmetry models like fli] , |45|j for instance, generically have a much larger 
Mi/M 2 hierarchy than we found necessary for leptogenesis. 

Second, one may question our crude implementation of the SO (10) relations 
between quarks and leptons Dirac- masses. We simply used mi ep ton = rriquark/3, which 
fails by about a factor 3 for the muon-strange mass ratio. Since the asymmetry 
goes like r 2 (see eq. |4.7|) , such small factors should not qualitatively change our 
conclusions; taking 10 instead of lMeV for the lowest neutrino Dirac mass could for 
instance barely bring the vacuum oscillation solution up to the required asymmetry 
Ybw ~ 1- The correct running of neutrino Yukawa couplings may have important 
effects on the mixing, which could change the interesting value of s e 3, and for instance 
make it inaccessible to first generation long baseline experiments, or on the contrary, 
already excluded by CHOOZ. 

Our conclusions seemingly contradict previous literature on the subject, which 
usually left the impression that leptogenesis could work without restrictions. This 
comes from our maybe over-restrictive framework, insisting on both 5*0(10) relations, 
and recent neutrino oscillations data. Let us review some of these discrepancies. In 
reference [ IS] , the authors use S'O(10)-inspired see-saw relations like eq. |3. 14 
and conclude that leptogenesis is generically possible. However, the values taken 
in the final analysis are hard to reconcile with maximal atmospheric mixing, which 
at the time was not as firmly established as today. In subsequent works [47|, f48[ 



maximal atmospheric mixing is well taken into account, but the Dirac mass pattern 
considered is derived from a family symmetry in an SU(5) framework; this gives 
typically larger "CKM" mixing and much weaker Dirac hierarchy than the simple 
quark-lepton ansatz we considered, both of which increase the asymmetry. In refer- 
ence it is shown that leptogenesis can be reconciled with SMA solar neutrinos. 
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But the Dirac masses needed badly violate 50(10) relations. Finally, another way to 
loosen the strong 50(10) constraints is to invoke what could be coined "see-saw in 
the scalar sector", namely a small v.e.v. for a scalar triplet that directly contributes 



to the left neutrino Majorana mass [31]. This however introduces a new parameter, 
and non-trivial constraints on the scalar sector as in |50 |. 



Note added: While this work was being revised, others |51|,|[52|| have confirmed 
the important effect of U e s in the framework discussed here. It should however 
be stressed that the type of random scan used in these works requires a very fine 
graining to hit the maximal possible value: at least 10 5 points in each f/ e3 plane, 
as seen figure |6|. After the recent KamLAND results [41], the only surviving LMA 
solution requires an optimal value of \U e s\ ~ 0.12 which is but slightly lower than 
the VAC case detailed in this work. However, reconciling thermal leptogenesis with 
an 5O(10) mass pattern seems even harder than before. 
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